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Abstract 

We study the dependence of quark condensate £ on an external magnetic field. 
For weak fields, it rises linearly: 



£(#) = E(0) 



eH]n2 „ f e 2 H r 
1 + ^ n + O 



16tt 2 F 2 



Fj 



(1) 



M n and F w are also shifted so that the Gell-Mann - Oakes - Renner relation is 
satisfied. 

In the strong field region, T,(H) oc (eH) 3 / 2 . 

1 Introduction. 



Phase structure of QCD is now a subject of an intense discussion. The question of how 
the properties of the system are modified by non-zero temperature was studied especially 
well (see e.g. [|J for a recent review). We know now that, in the theory with massless 
quarks, phase transition with restoration of spontaneously broken chiral symmetry occurs 
at some temperature T = T c . That means that the order parameter of the symmetry 
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breaking, the quark condensate £ = — < qq > falls down as temperature increases and 
turns to zero at T > T c . When temperature is small compared to characteristic hadron 
scale fihadr, the dependence S(T) is known exactly |2|. In the case of two massless flavors, 



S(T) = E(0) 1 



8F2 384F4 



(2) 



The derivation of this formula relies on the fact that a lukewarm heat bath involves 
mainly pions — other degrees of freedom are not excited yet. The pion interaction at 
small energies is known from the effective chiral Lagrangian: 



Here U is a unitary SU(2) matrix ( it may be parameterized as U — exp{ir a a /F T } 
where <j) a is the pion field) and we included also the mass term involving quark mass matrix 
M.. F v = 93 MeV is the pion decay constant and the parameter £ has the meaning of 
quark condensate H = \ < uu > \ = \ < dd > \. The particular formula (|j) of the 
Lagrangian is dictated by chiral symmetry (see || for a nice pedagogical review). 

But temperature is not the only external parameter which can affect the properties of 
the system. One can consider equally well a cold but dense system with non-zero mean 
baryon charge density related to the chemical potential. This system is less studied, but 
there are good reasons to believe that also in this case chiral symmetry is restored when 
the chemical potential exceeds some critical value. 

In this paper, we address the issue of how the properties of QCD vacuum state de- 
pend on an external magnetic field. A naive expectation based on the analogy with 
superconductivity and with the situation in QCD at non-zero temperature and/or chem- 
ical potential could be that the condensate decreases as the magnetic field increases and 
melts down completely at some critical value H = H c above which the chiral symmetry is 
restored. We will see (and that is our main conclusion) that it is not so. Quark condensate 
rises with the increase of magnetic field and no phase transition occurs. 

The behavior of a hadron system in magnetic field was studied earlier by Klevansky 
and Lemmer in the framework of Nambu-Jona-Lasinio (NJL) model [Q. Solving gap 
equation in the presence of an external magnetic field they found that the order param- 
eters of the spontaneously broken chiral symmetry, the dynamical fermion mass and the 
chiral condensate, rise with the field. This conclusion was confirmed in later studies of 
NJL model and related theories ||. The shift of fermion mass and of the condensate is 
quadratic in field 



So, in this model, chiral symmetry is not restored by magnetic field. Qualitatively, 
NJL method behaves in the same way as QCD. However, we shall see later that the 
quantitative predictions are different. 



L = -Z-TkidfJJdvtf} + SRe Ti{MU ] } + higher order terms 



(3) 
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Figure 1: Vacuum energy in a weak magnetic field. Solid lines stand for charged pions. 



2 Weak field. 

We will consider QCD with two massless flavors. In the leading order of chiral perturba- 
tion theory, the masses of u— and d— quarks appear in the vacuum energy only via its 
dependence on M% which is proportional to the sum m u + m d - It means that in this order 
< uu > and < dd > condensates will not differ and we can define the quark condensate 
as 



wm de vac (m u ,m d ,H) 
H H ) = 



dm, 



vac 

(m u ,m d ,H) 

m u =m l {=0 



dm d 



(5) 

=m d =0 



where the small quark masses m u , m d are introduced as external probes. 

To find e vac , we need to calculate vacuum loops in the presence of magnetic field. When 
the field is weak eH <C fi 2 adr ~ (2irF n ) 2 (however, we will always assume eH ^> M^ 2 ), 
characteristic momenta in the loops are small and the theory is adequately described by 
the effective low energy chiral Lagrangian (|3]). In the leading order, pion interactions can 
be neglected whatsoever, and the field dependent part in e vac is given by one loop graphs 
depicted in Fig. 1. 

It is instructive to consider first the graph with two photon legs. The corresponding 
contribution in the vacuum energy is 

= (6) 

where A is the ultraviolet cutoff. By virtue of Gell-Mann - Oakes - Renner relation, 

F 2 M 2 = Z(m u + m d ), (7) 

one can relate the derivatives over m u {m d ) and over M n 2 . Taking into account Eqs. (^), 
(H), and (0), we obtain 

A^ 2 \H) = — 2 6 f 2 - (8) 

V ; F n 2 96n 2 M n 2 V ; 

This expression diverges in the chiral limit M w 2 — > and makes as such little sense. 
It is easy to trace back the origin of this divergence. The corresponding graph for e vac 
involves the logarithmic divergence both in the ultraviolet and in the infrared. After 
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differentiating, it gives the contribution in the condensate which diverges as a power in 
infrared. 

But one is not allowed, of course, to restrict oneself by the graphs with only two 
external field insertions. All other graphs are also important. The more is the number of 
legs, the more severe are the infrared singularities. For example, the graph with four legs 
involves the sing ularity ~ (eH) A /M n 4 in the vacuum energy which gives the singularity 
oc l/M n 6 in the condensate, etc. All such graphs should be summed up. The result is the 
analog of the Euler-Heisenberg Lagrangian for scalar particles which was found long time 
ago by Schwinger ||. We have [] 



167T 2 jo s 



00 ds 



-Mtt 



eHs 



smh(eHs) 



- 1 



(9) 



We see that the integral (H) is regular in the infrared (large s region). The infrared cutoff 
is provided now not by M^ 2 , but by the field eH itself. All infrared divergent pieces in 
the graphs in Fig. 1 are summed up into a finite expression. 

To find the shift of the condensate, we have to differentiate Eq. @ over quark mass 
and substitute it into the definition (||). We arrive at the expression proportional to the 
simple integral 

z l [ smh{z) J 

which can be done as a sum of residues of the poles on, say, positive imaginary z axis. 
Our final result is 

E(/ ^ s(0) [ 1+ ii!^ + 'i <»> 

We see that the shift is positive and linear in H. The latter is easy to understand if 
substituting the actual infrared cutoff ~ eH for M 2 in Eq. (|8|). 

It is also clear now why NJL model gave the shift which was quadratic rather than 
linear in field: the corresponding calculation involved the loop of massive quarks rather 
than the loop of massless pions and was infrared finite in any order in H. A side remark 
is that a simple-minded NJL calculation is also not able to reproduce the result (|2|) for 
the temperature dependence of the condensate. At small temperatures, the density of 
massive quarks in the heat bath is exponentially suppressed, while massless pions are in 
abundance. 

The correction cx eH in Eq. (|TTJ) is of order 1 when the field is comparatively large 
V eH ~ 1.4 GeV. Well before that, our calculation based on the effective pion Lagrangian 
(ID loses validity. In principle, one can trace the deviation from the leading order result 
Tl"D at intermediate values of H in the framework of chiral perturbation theory p|] taking 



into account nonlinear pion interactions. An example of the graph contributing in the 



lr rhis is the unrenormalized vacuum energy. More precisely, the quartic divergence which corresponds 
to the loop without legs and is field independent is subtracted, but the logarithmic ultraviolet divergence 
associated with the graph with two legs is not. The expression usually found in the textbooks corresponds 
to subtracting also the two- leg graph, so that the expansion of e^(H, M„) in H starts from the term 
oc H 4 . The two-leg graph is then absorbed in the charge ( = field ) renormalization, and the derivative 
e^|£(iJ, M„) over mass would, of course, be the same as that of 



of H? en (M„)/2 



c unrcn 

-vac 
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Figure 2: A two-loop graph contributing to the vacuum energy. 

vacuum energy in the next order is depicted in Fig. 2. Also higher order terms in the chiral 
Lagrangian which we did not specify and did not discuss here would become important. 
The estimate for two-loop correction to the condensate is oc £(0)(eif ) 2 / (27rF n ) 4: (the same 
estimate oc S(0)(e/f) 2 / / u^ adj , is obtained if taking into account massive charged particles 
such as K and p meson). In this order, one should not expect the corrections to | < uu > | 
and | < dd > | to be the same (the magnetic field breaks down isotopic invariance). The 
calculation of the coefficients (it is not yet clear whether also some logarithmic factor 
oc In appears ) is now in progress. 



3 M W (H), F W (H) and Gell-Mann - Oakes - Renner re- 
lation. 

As the electric charges of u— and d— quarks are different, flavor symmetry is broken in an 
external magnetic field. In particular, the axial SUa{2) symmetry is broken down to U A (1) 
corresponding to chiral rotation of u— and d— quarks with opposite phases ( the singlet 
axial symmetry is broken already in the absence of the field due to the anomaly ). The 
formation of the condensate breaks down this remnant U\{1) symmetry spontaneously 
leading to appearance of a Goldstone boson, the 7r°-meson. Indeed, charged pions acquire 
a gap in the spectrum oc \J eH and are not goldstones anymore. If quarks are endowed 
a small non-zero mass m, pions are not exactly massless, their mass being related to the 
quark condensate by the Gell-Mann - Oakes - Renner relation ([?]). 

It is interesting to study the question how the mass and residue of 7r° depend on the 
external field. 

Let us find first the mass shift. We have to calculate the polarization operator of 7r° 
given by the sum of graphs of the kind drawn in Fig. 3. 

The four-pion vertex can be found from the expansion of the effective chiral La- 
grangian. In the exponential parameterization U = exp{i<fr a T a / F n }, the corresponding 
terms in the Lagrangian have the form 

L {i) = ^[(<W a ) 2 - (W(<9 M b ) 2 ] + ^(W 2 - (12) 

We need also the expression of charged pion propagator in a magnetic field. It can be 
inferred from the results of Ref. || where an explicit expression for the fermion propagator 
in a magnetic field at non-zero chemical potential \x has been found. Basically, one has to 
take the integral multiplying the factor — //7 in Eq. (4.9) of Ref. || without the factor 
cos(eifs) + 7 X 7 2 sin(eifs) = exp{ie(aF)s/2} in the integrand. We arrive at the following 
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Figure 3: An example of the graph contributing to the polarization operator of 7r° in a 
magnetic field. 



expression for the Euclidean scalar propagator 

D B (x,y) = exp{ie J* A M / ^e~ lk ^D H (k), 



d A k 



(27T) 4 



(13) 



where the integral in the phase factor is done along the straight line connecting x and y, 
and 



D H {k) 



ds 



■ exp 



k 2 + k'' 



tanh(eiJs) 
- eHs 



+ M1 



(14) 



'o cosh(eHs) 

with k?, — kl + k^, kj_ — kf + k\ (the magnetic field is aligned along the third axis). 

Substituting (|i~4"|) and the vertex inferred from ([12]) in the graph in Fig. 3 (as the 
propagator enters with coinciding initial and final points, the phase factor disappears), 
and subtracting the similar expression at H = 0, we obtain 



iW)-iM0) = ~f -^(A4 2 +2pW) r ds 

3F* J (2-7T) •'0 



exp{- S (^ + A;i^M£) + M 2 )} 



cosh(eifs) 



(15) 



The mass shift is given by the shift fllq ) of the (Euclidean) polarization operator at 
the point p 2 = —M 2 . Note first of all that the shift is zero for massless pions. An exact 
goldstone remains the exact goldstone also when H ^ 0. The shift A^M 2 is proportional 
to M 2 (H = 0). Taking two first terms of the expansion of M 2 and calculating first the 
integral over d A k = ir 2 dk 2 dk\ and then over ds (the latter has the same structure as 
the integral ( ]TD| ) which we encountered when calculating the shift of the condensate and 
involves also an extra piece ~ J °° ds[cosh(s) / sinh 2 (s) — 1/s 2 ] which is, however, zero), we 
arrive at the simple result 



„ , 9n „ , 9 eifln2 „,r,ei?ln2 



247r 2 F 2 T T 487r 2 F 2 



+ 



which implies 



Ml{H) = M 2 (0) 



eH In 2 
1 ~ 16n 2 F3 + 



(16) 



(17) 
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Figure 4: Axial currents correlator in a magnetic field. The graphs contain exact propa- 
gator of charged pions in the constant magnetic field. 
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Let us find now the renormalization of the residue brought about by the field. We need 
to calculate the one-loop graphs depicted in Fig. 4 which contribute to the pole structure 
in the axial current correlator / dxe ipx < A 3 (x)Al(0) > H d 4 x ~ p fl p u F^(H)/(p 2 + M%). 
As F n and its shift are non-zero in the chiral limit, we can set M 2 = 0. The graph in Fig. 
4a has already been calculated (the first term in Eq. fliTf ) ). To calculate the graphs in 
Fig. 4b, c , we need the vertex < 0|v4^|7r 7r + 7r~ > which can be found by " covariantizing" 



the derivative in Eq. (^) according to the rule [0] dyU — > d^U — ^A^U + UA^) where 



A^ = A^t a . The final result is 



eH\n2 
1 + 8tt 2 F 2 + 



Note that the contributions of the individual graphs Fig. 4a and Fig. 4b, c depend on 
the chosen exponential parameterization, but the physical pion residue F n (H) depend- 
ing on their sum does not. Also the physical pole position (|I7D does not depend on 
parametrization while the individual terms in flTE| ) do. 

The results (p|) , ([l7|) , and fllBD are very much analogous to Eq. (||) and the similar 
known expressions for M%(T) and F%(T). Only here the physical situation is just opposite 
compared to the thermal case. A nonzero temperature suppresses the condensate and 
residue, and leads to the increase of the mass. A nonzero magnetic field brings about the 
increase in the condensate and residue, and suppresses the mass. 

In the thermal case, the renormalized condensate, pion mass, and residue still satisfy 
the relation (0) in the leading order in T 2 . Likewise, the expressions ([LTD >(H3)j an d flI5|) 
satisfy the Gell-Mann - Oakes - Renner relation in the leading order in eH. 



4 Strong field. 

When eH 3> ^\ a dri characteristic momenta in the vacuum loops are high, and the system 
is adequately described in terms of quarks and gluons rather than in terms of pions and 
other low lying hadron states. Due to asymptotic freedom, the effective coupling constant 
a s (eH) is small, and we can try to treat strong interaction effects pertubatively. 

For sure, in QCD with massless quarks, the condensate cannot appear in any finite 
order of perturbation theory, irrespectively of whether a magnetic field is present or not 
- both electromagnetic and strong interaction vertices respect chirality. 

It was recently discovered, however, that the condensate is generated in the strong field 
limit ||. To see that, one has to sum up an infinite set of relevant graphs. The authors 
of || studied the Bethe-Salpeter equation which implements such a resummation for 
massless QED and showed that the equation admits a nontrivial solution with dynamically 
generated mass. Later this result was reproduced in the language of Schwinger- Dyson 



equation ||10|| . The derivation applies also in QCD without essential modifications. For 
clarity sake, let us say here few words about it. 

When strong interaction is disregarded whatsoever, we have a system of free charged 
quarks in a magnetic field. The spectrum presents the set of Landau levels 



e ± (n,a,k 3 ) = ±J\e q H\(2n + a + l) + kl + ml (19) 
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where n = 0, 1, a = ±1 marks the spin orientation, and k 3 is the momentum along 
the magnetic field direction. Negative energies describe the Dirac sea. The corresponding 
eigenstates are localized in the transverse direction. The characteric size of the orbits is 
~ \j\feH ( so that for \JeH ~ 1.4 GeV they are already pretty small). In infinite space, 
each state (|IP|) is infinitely degenerate, the degeneracy being associated with the position 
of the center of the orbit. For a finite box of size L, the level of degeneracy is 

N ± = ^L 2 , (20) 
Ztt 

and also k 3 is quantized to 27rn/L. For m q = 0, the ground states have zero energy. 
We need in the following the fermion Green's function in a magnetic field. Its ex- 



plicit expression was found in [11] and, in a somewhat more convenient form, in M with 



Schwinger technique. First of all, one can write 

G H (x,y) = exp |ie 9 J A^Qd^ G H (x - y) (21) 

where the integral in the phase factor is done along the straight line connecting x and 
y. The Fourier image of Gh{x — y) presents a complicated integral over proper time. 
Fortunately, we do not need the full expression, but only its asymptotic form in the 
region of small momenta k <C \feH. In this region, it suffices to retain only the lowest 
landau levels (LLL) with n = 0, a = — 1 in the spectral decomposition of the Green's 
function, and the latter acquires the simple form [11], M 



G H {k) = ie H*/l««*l> 1^^1(1 - ryY), (22) 



where fcjj = k% — k\ and k\ = k\ + k\. 

Basically, this Green's function describes a free motion of the states with a = — 1 in 
longitudinal direction. Strictly speaking, retaining the exponential factor exp(— k\/\e q H\) 
is not quite consistent - the dominance of LLL on which the derivation of Eq. ( p2|) was 
based is justified only in the region where both k\ and k^ are small compared to \e q H\ 
and the exponential factor is not efficient. We will see later, however, that the generation 



of condensate is related to the infrared region k J\e q H\. Really, in the opposite limit 



k ^> J\e q H\, the Green's function Gh(x, y) tends to the free fermion Green's function, 
and we cannot expect a nontrivial dynamic phenomenon like the condensate generation 
to be associated with that region. The exponential factor in (|22|) is convenient to retain 
as an effective momentum cutoff. 

We are interested in the theory with massless quarks, so that m q = in the first place. 
We expect, however, the dynamical mass generation, and our Ansatz for the exact Green's 
function at low momenta is Eq. ( P^ ) with nonzero and, generally speaking, momentum- 
dependent m q . 

Next, we substitute this Ansatz in the Schwinger- Dyson equation schematically pre- 
sented in the Fig. 5. The equation was solved in the approximation where the loop 
corrections to the vertices and also to the gluon propagator (the latter is actually a rather 
strong and not so innocent assumption) were disregarded. 
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Figure 5: Schwinger-Dyson equation for the exact quark propagator. Thin solid lines 
correspond to the free Green's function and bold lines - to the exact one. 



Doing the spinor trace and performing Wick rotation, one can see after some trans- 
formations that the equation admits a factorized solution 

™> q (p\\,p±) =rhq(P\\ 2 )exp{-P±/\e q H\}, (23) 
where m 9 (p|| 2 ) satisfies the following 2-dimensional integral equation 

A 2 a s c F f d 2 k\\ m q (k\\ 2 ) f d 2 k ± exp{-kj/(2\e q H\)} 
m ^ ) = ^J kl + ra^) J (kn-puf + fci • (24) 

The integral (Q) is double logarithmic f\ receiving a main support from the region 
m 9 2 (0) < k 2 < k]_ < \e q H\ 

Let us first tentatively neglect the dependence of rh q on longitudinal momenta. We 
would obtain 

m g (0)~m 9 (0)^ln 2 ^§. (25) 
47r m 2 (0) 

The corresponding calculation was actually first done in [^] where a correction to the 
electron mass in the presence of the external mag netic field Am (H) = ^ ln 2 (e#/m 2 (0)) 
was found. In that (and other) earlier papers, the equation (^) was not treated, however, 
as a self-consistent equation allowing to unravel the dynamical generation of mass even if 
the lagrangian electron mass is zero. The solution of this equation is 



SW-v^expj-^] (26) 

The exponential factor displays a truly nonpertubative nature of the result. An accurate 
analysis of Ref. || which takes into account the momentum dependence of the mass leads 
to the result 



7T 7T 



A ~ V |e,i/|e X p|--y— | (27) 
The quark condensate is defined as 

< qq > H = - / ——4 TrG H {p) = -77—4 / d P± d P\\ 2 , 2 , 2( 2 2 \ ( 28 ) 
J (27r) (2-k) J P± +pj + ra 2 Gojj,Pi) 



2 The double logarithmic structure of the integral is easy to understand if disregarding the phase 
factor in Eq. (|2l| ) which makes the calculations elementary. That gives qualitatively the same form of 
the integral as in Eq. (|24|). Quantitative results, however, would be wrong. In particular, the coefficient 
in Eq. (|2^) would be 4 times larger than the actual one. 
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Assuming the exponential fall-off of dynamical mass at p\ ^> \e q H\, p| 3> \e q H\ (cf. Eq. 
fl23D ), we finally obtain 

S H = - < ft >„ = ^^.p ^^ l (29) 

where F(a s ) is some function not involving an exponential dependence. 

We see that the condensate increases with the field, and no phase transition occurs. 



5 Discussion. 



Our main result ( |TT]) is an exact theorem of QCD. It has the same status as Eq. (|2]) and 
some exact results for density of eigenvalues of Euclidean Dirac operator at finite volume 
TM and in thermodynamic limit |[L5||. It would be very interesting to check this and other 



similar exact theorems in lattice calculations with dynamical fermions and/or vacuum 
model calculations Q. 



The derivation of the result (|29 ) in the strong field limit reviewed in the previous 
section is similar in spirit to the derivation of spontaneous symmetry breaking in the NJL 
model. In both self-consistent gap equation for the fermion propagator is solved. 

This method of derivation seems to be rather reliable , but it is indirect. It would be 
very nice to understand the appearance of the quark condensate in strong fields in plain 
and direct terms. 

We have in mind the following. According to the famous Banks and Casher theorem 



JT7[ , the quark condensate is related to the spectral density of Euclidean massless Dirac 
operator p(A) at A = 

£ = 7rp(0). (30) 

Nonzero spectral density p(0) means that the number of Dirac operator eigenvalues with 
A < Ao fihadr is estimated as iV(A < Ao) = p(0)XoV, where V is the Euclidean volume. 
For free massless fermions, the spectrum is A n = 27r|n|/L where n is a four- dimensional 
integer vector. p/ ree (A) ~ A 3 , p/ ree (0) = and the condensate is zero. For free massless 
fermions in a magnetic field, the infrared branch of the spectrum is 

2tt, . , , 

A m = — |m|, (31) 

where m is a two-dimensional integer vector. The eigenvalues (|3T|) correspond to free 



motion in the longitudinal direction. Each level (pl|) involves the high level of degeneracy 



(|). The spectral density is then estimated as pf ree (X) ~ \e q H\\. This is "better" 



3 Unfortunately, lattice calculations are very difficult here. The Euclidean lattice should be roomy 
enough to accommodate pions: the size of the box should be considerably larger than the pion Compton 
wavelength. Calculations in the framework of the instanton liquid model ||l6f allow for larger boxes and 
are more promising. 

4 There is, however, a not quite clear for us at the moment question on whether the loop corrections to 
the gauge boson line in the Schwinger - Dyson equation in Fig. 5 could modify the result. According to 
|l3} | , the effective mass which gauge bosons acquire in a magnetic field brings about an additional infrared 
cutoff which kills one of the logarithms in the strong field asymptotic of the fermion mass operator 
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than in the absence of the field, but still E = 7rp(0) = 0. A nonzero condensate ( p9|) 
means that, after taking into account the interaction between quarks, the Euclidean Dirac 
spectrum (|3TD is substantially modified, and small eigenvalues with characteristic spacing 



~ \e q H\L 2 ) appear. We are not able now to display the presence of such small 
eigenvalues explicitly. 

An external magnetic field increases the condensate which means that it should make 
the chiral restoration phase transition in temperature and/or in baryon chemical potential 
more difficult. That means, in particular, that critical temperature T c (at H = 0, it is 
estimated to be of order 200 MeV |TJ) should increase with H. This question was studied 



in recent [18] for massless QED. Translating their result in QCD language, we obtain 



T c ~a s J\e q H\ (32) 



in the strong field region. The result fl32|) is easy to understand. J\e q H\ appears by 
dimensional reasons and the factor a s appears because the generation of the condensate 
in a strong magnetic field is driven by strong perturbative inter-quark interaction. 
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